Singlet-triplet relaxation induced by confined phonons in nanowire-based quantum 
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The singlet-triplet relaxation in nanowire-based quantum dots induced by confined phonons is 
investigated theoretically. Due to the quasi-one-dimensional nature of the confined phonons, the 
singlet-triplet relaxation rates exhibit multi-peaks as function of magnetic field and the relaxation 
rate between the singlet and the spin up triplet state is found to be enhanced at the vicinity of the 
singlet-triplet anti-crossing. We compare the effect of the deformation-potential coupling and the 
piezoelectric coupling and find that the deformation-potential coupling dominates the relaxation 
rates in most cases. 



I. INTRODUCTION 

As one of the most promising candidates for qubits, 
spins in semiconductor quantum dots (QDs) have at- 
tracted much attention in the last two decades*^— Al- 
though much progress has been made in single-electron 
QDs, an attractive alternative is to base the qubits 
on singlet-triplet (ST) states in two-electron QDs4^— 
This is mainly due to two important features of the ST 
states. One is the Zeeman-driven ST transition with anti- 
crossing due to the spin-orbit coupling (SOC), which en- 
ables coherent manipulation of spin states^— The other 
one is the low ST relaxation rate which makes the spin 
polarization remains for a sufficiently long timeJ^— As 
the QDs are either self-assembled ones sitting on the sur- 
face of the substrate or fabricated by confining electrons 
in quantum wells by electrodes, bulk phonons in the 
substrate play an important role. In low-temperature 
regime, acoustic bulk phonons in conjunction with the 
SOC serve as the main source of the ST relaxation due 
to the suppression of the effect of the hypcrfine inter- 
action in the presence of large magnetic field. Previous 
studies show that at the vicinity of the ST anti-crossing, 
the relaxation rate is greatly suppressed due to the mis- 
match of the phonon emission wave length to the dot size. 
Furthermore, the spin-up triplet states which couple to 
the singlet ground state through the SOC has a much 
shorter lifetime compared to the other two triplet states 
due to the strong ST mixing between themJ^r— 

Despite of these progress in two-electron QDs, the com- 
plicated fabrication process makes them difficult to be 
scaled up, which limits their application from an in- 
dustrial point of view. Recently, self-assembled two- 
electron QDs embedded in InAs [111] nanowires were 
fabricated^ - — The well-controlled growth process en- 
ables parallel production in massive number with simi- 
lar properties, which makes these nanowire-based QDs 
are more suitable for integration onto silicon and scal- 
ing up for large hierarchical systems. Although simi- 
lar Zeeman-driven ST transition has been obtained in 
nanowire-based two-electron InAs QDs^ir— the ST re- 



laxation can be quite different due to the properties of the 
confined phonons in the nanowires. The nanowires are 
perpendicular to the substrate, making the bulk phonons 
in the substrate are less important than the confined 
phonons in nanowiresi 27 ' 28 The regular structure of the 
nanowires results in the quasi-one-dimensional confined 
phonons, which lead to novel properties in optical ab- 
sorption and transport for the nanowire-based QDs i 29 i 30 
Therefore, these confined phonons are expected to have 
a pronounced impact on the ST relaxation, which has 
not been studied to the best of our knowledge. In this 
work we will investigate the ST relaxation induced by the 
confined phonons. 



II. MODEL AND FORMALISM 

Consider a two-electron elongate QD embedded in 
InAs [111] cylindrical nanowirc with radius R in the 
presence of external magnetic field B along the wire. 
We model the QD by an anisotropic harmonic poten- 
tial V c (r,z) = ^m*u> 2 r 2 + \vh*uj 2 z 2 ( corresponding to 



the effective dot length d z = y/tm /m*(J z and dot diam- 
eter do = -\/?nr/m*ujo) with z-axis along the wire, m* 
is the effective mass of the electron. The single-electron 
Hamiltonian can be expressed as 



H c = 



_E_ 

2m' 



+ V c (r,z) + H B + H SCh 



(1) 



where H& = ^g/j,^B ■ cr is the Zeeman splitting with g, 
/Ub and cr being the g factor of electron, Bohr magne- 
ton and Pauli matrix respectively, -ffso represents the 
SOC term. In this work, we concentrate on the Rashba 
coupling which is dominant in InAs nanowires which has 
the form if go = \& V V Z with 7 being the Rashba cou- 
pling strength.— i^i We assume the dot diameter do <C d z 
so that only the lowest electron subband in the radial 
direction is needed. We also assume the external mag- 
netic field is weak enough so that its orbital effect can be 
neglected. 



The confined phonons arc calculated with isotropic 
elastic continuum model which is widely used in the study 
of nanowires, carbon nanotubes and nanoparticles.— ~— 
The nanowire is modeled as an infinite cylinder with the 
stress vanishing at the surface of the wire. The displace- 
ment field of the ions induced by the confined phonons 
can be expressed by the quantized form (in cylindrical 
coordinate) 



u(r) 




u vq {r)e u °^{a vq + al_ q ), (2) 



where V is the volume of the nanowire and p is the den- 
sity. The eigenmodes u ua (r) and spectrum tu uq are calcu- 
lated following Rcfs. I38ll39l . Note the due to the bound- 
ary condition, generally speaking, the confined phonon 
modes have both longitude and transverse components. 
Given the displacement field, both the deformation- 
potential coupling and the piezoelectric coupling can be 
calculated. The deformation-potential coupling is given 
by the divergence of the displacement field, which reads 



H»Jr) = -ZV-u(r), 



(3) 



where 5 is the deformation-potential coupling strength. 
Note that only the longitude component of the con- 
fined phonon mode has contribution to the deformation- 
potential coupling. The piezoelectric coupling is given 
by^ 



H?Jr) 



dr. 



V-P 



PZ 



(4) 



where P pz is the polarization induced by the displace- 
ment field and k is the dielectric constant. For InAs [111] 
nanowires with wurtzite structure, P pz can be expressed 
asia 

P pz (r) = e 15 (d r u z (r) + d z u r (r)\e r 

+ e 31 (d r u r (r) + u r (r)/r) + e 33 d z u z (r) e z , (5) 

with ei5, e 3 i and 633 being the piezoelectric constants. 
Note that due to the anisotropy of the polarization, the 
piezoelectric coupling can be more sensitive to the con- 
fined phonon mode than the deformation-potential cou- 
pling. 

In the following discussion, we restrict to the dilatation 
modes [which have angular momentum quantum number 
I = and angular component u B v (r) = in Eq. ([2"jl] since 
only these modes can couple to the electrons in the QDs 
in the lowest subband of radial directioni 30 ' 41 i 42 

For two-electron QDs, the total Hamiltonian of the 
system is given by 



H = (Hi + H 2 + H C ) + Hl p + H 2 cp + H p , 



(6) 



where Hq 



ki-r 2 | 



represents the Coulomb interac- 



tion between the two electrons. H p = VJ hxi vq a[, a vq 



is the Hamiltonian for the confined phonons and H op = 
H® p + H P p is the electron-phonon coupling. The super- 
script "1" and "2" label the two electrons. 

To construct the two-electron basis functions, it is con- 
venient to use the separation of variables in terms of the 
center of mass R = (ri + r2)/2 and the relative motion 
r = t*i — T2- As the Coulomb interaction is too strong, 
we choose the basis function as the eigenstates of the 
Hamiltonian (H\ + H 2 + He) without SOC and Zceman 
term, which can be solved numerically by finite difference 
method in real space. The basis functions have the form 

(R,r\n R n, V ) = ^ R (i?)^ r (r)^ R (Z)^„ r (z) Xv , (7) 

where \n represents the spin states of the two electrons, 
which can be expressed as 



X V 



1 5 ) = 75(1 n> - 1 n», v = q 

|T_> = |U>, V = l 

|T > = ^(in) + in», v = 2 

|T+) = |tt>, 77 = 3 



(8) 



In Eq. (0, 7Z R (R)/7Z r (r) is the wave function for the low- 
est subband in radial direction and <$% R (Z) / 4> r (z) is the 
wave function in axial direction. The subscript "iZ" and 
"r" label the center of mass and relative motion respec- 
tively. To guarantee the anti-symmetric of the wave func- 
tion, the singlet state \S) always corresponds to <f? n (z) 
with even parity, while the triplet states iTi^) corre- 
spond to <p T n (z) with odd parity. The energy levels and 
eigenstates of the two electrons can be obtained by diago- 
nalizing the full two-electron Hamiltonian (H^+H^+Hc) 
in these basis. We identify an eigenstate as singlet and/or 
triplet by its expectation value {(<t\ + <r 2 ) 2 ). 

Treating \i) and |/) as the initial and final states, we 
can calculate the phonon-induced relaxation rate using 
Fermi's golden rule. At zero temperature, the relaxation 
rate induced by the confined phonons reads^ 



fi 



E 



\M 3 

I vqi 



Af\v2«n\*)\ 



\d qm (fi^u(.q m ))\ 



-0(E) 



(9) 



huj u (g„ 



where E = \Ef — Ei\ is the energy splitting and 
9(E) is the step function. Vj q = Wi q (r x )e lqzi + 
Wl (r2)e tqZ2 comes from the total electron-phonon in- 
teraction Hamiltonian H^ + H^ p , with j = D for 
the deformation-potential coupling and j = P for 
the piezoelectric coupling. The corresponding coeffi- 
cients are \M^ q \ 2 = KE? /(2t< pu vq R 2 ) and \M* q \ 2 = 

\Qh-K 2 e 2 e\ A /(n 2 2-KR 2 uj vq ). The quantity W uq (r) are 
given in detail in Appendix. Note that the relaxation rate 
induced by bulk phonons can be expressed with similar 
equation as Eq. ([§]), the main difference is that since the 
constant-energy surface is continuous for bulk phonons, 
the discrete summation over q m should be replaced by 
the integration over wave vector q. 

Before we represent the results, we would like to briefly 
summarize the three mechanisms which is crucial for the 



ST relaxation. The first one is the ST mixing. The ST re- 
laxation favors strong ST mixing where the SOC couples 
the states more efficiently. The second one is the ratio be- 
tween the phonon emission wave length and the QD size. 
The phonon emission efficiency is maximized as the two 
length is comparable to each other, leading to large relax- 
ation rate. The third one is the phonon density of states 
(DOS) at the phonon emission energy, which is equal to 
the ST energy splitting. Larger DOS implies that more 
phonons can be emitted, which enhances the ST relax- 
ation. The first two mechanisms manifest themselves in 
(/ \Vj | i) in Eq. ©, which is known as the form factor. 
The third one, e.g., the phonon DOS, is closely related 
to the quantity d qm (ftiu v (q m )), which is very sensitive to 
the phonon spectrum. Thus this one is expected to lead 
to pronounced differences in the ST relaxation induced 
by confined phonons and bulk phonons. 



10 



ID — - - 



W 9 - 



n - 





~~X- 


_ 


X 


V 


~ 



12 3 4 

B(T) 

FIG. 1: (Color online) The lowest 4 energy levels and expec- 
tation value of S~ = (<rf + o"i)/2 v.s. external magnetic field 
B. 



III. RESULTS 

In the numerical calculation, we use the parameters for 
the InAs [111] wurtzite nanowires i ~ 26 The deformation- 
potential strength S is chosen to be 5.8 cV, the density 
is p = 5900 kg/m 3 , the static dielectric constant k is 
15.15, the longitude and transverse sound velocities are 
chosen to be vl = 4410 m/s and vt = 2130 m/s re- 
spectively^ Since the piezoelectric constant for wurtzite 
InAs nanowire is largely missing, we use the constants 
calculated by transformation from zinc-blende value en, 
which are eis = e^i = —633/2 = —en/v3, with 
eu = 3.5 x 10 8 V/mM& The g factor is set to -9.0,21 
and the Rashba coupling constant 7 is chosen to be 
1.5 x 10~ n eV-m (corresponding to the SO length Aso = 
200 nm)4£ We set the radius of the nanowire R to 15 nm. 
The effective QD diameter do is set to 8 nm. The dot 
length d z is set to 70 nm. 

We employ the exact diagonalization method with the 
lowest 48 basis functions to converge the energy levels 
and ST relaxation rates i " The magnetic field depen- 
dence of the first four levels and the corresponding expec- 
tation value of S z = (a\ + erf )/2 are plotted in Fig. QTa) 
and (b) respectively. An anti-crossing at B = 2.36 T can 
be identified in the energy levels. The expectation val- 
ues (S z ) show an crossing at the corresponding magnetic 
field, indicating a strong mixing between the singlet and 
spin up triplet state (ST mixing). As the orbital effect 
of the magnetic field is neglected, the levels show linear 
field dependence, which is different from the typical disk- 
shaped QDsi 24 ' 31 We emphasis that the mixing between 
the singlet and triplet states in our elongate QDs is more 
stronger than the QDs studied before i 15 ' 17 ' 20 This can be 
seen in the large level splitting, which is 0.14 meV at the 
anti-crossing. This value is much larger than the typical 
value for disk-shaped GaAs QDs, which is usually a few 
/ieV. The strong mixing is not only due to the strong 
SOC in InAs, but also due to the large Coulomb interac- 
tion in the elongate QDs, which can enhance the relative 



strength of SOC H2L 

We start our discussion with the ST+ relaxation. The 
magnetic field dependence of the ST + relaxation is shown 
in Fig. [2ja) . It can be seen that the relaxation rates in- 
duced by bulk and confined phonons are very different 
from each other. For bulk phonons, the relaxation rate 
varies smoothly as a function of magnetic field. A local 
minimum can be found at the anti-crossing. For the con- 
fined phonons, the relaxation rate exhibits multi-peaks 
as a function of magnetic field. The two "sharp" peaks 
located at B = 1.11 T and B = 3.62 T are actually two 
divergent peaks. Moreover, the relaxation rate exhibits 
a local maximum at the anti-crossing in contrast to the 
local minimum for the bulk phonons. We attribute the 
multi-peak structure and the local maximum at the anti- 
crossing as the two main features of the ST + relaxation 
rate induced by confined phonons. 

The first feature, e.g., the multi-peak structure of the 
relaxation rate, is a direct consequence of the large DOS 
at the van Hove singularities in the confined phonon 
DOS. To see this, we plot the spectrum of the rele- 
vant confined phonon modes and the corresponding DOS 
in Fig. [2jb) and (c) respectively. The ST+ relaxation 
induced by the confined phonons via the deformation- 
potential coupling and the piezoelectric coupling arc plot- 
ted as a function of energy splitting in Fig. EJd) [Here 
we only plot the relaxation rate before the anti-crossing 
for clarification, the relaxation after the anti-crossing has 
similar behavior]. It is easy to see that the peaks in the 
relaxation rate correspond to the van Hove singularities 
in the DOS, indicating the strong enhancement due to 
the large DOS. The divergent peaks correspond to the 
van Hove singularity with q 7^ 0, while the non-divergent 
peaks are induced by the van Hove singularity with q = 0. 
This is because as q — > 0, the form factor tends to zero, 
which will suppress the relaxation rateil 3 - Note that the 
form factor can also have zeros for g ^ 0, which results 
in the dips in the relaxation rate. There are two types of 
zeros. One is the zeros due to the axial component (cor- 
responding to common peaks for both the deformation- 
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FIG. 2: (Color online) (a) ST+ relaxation rate as a function 
of external magnetic field. The black solid/red dash-dotted 
curve represents the relaxation induced by confined/bulk 
phonons. The thin yellow line indicates the position of the 
anti-crossing, (b) The energy spectrum of the relevant di- 
latation modes. Black solid/Blue dotted curves represent the 
axial/radial mode as q — > 0. (c) The corresponding phonon 
DOS. The arrow indicates the van Hove singularity with 
phonon wave vector q ^ 0. (d) ST+ relaxation rate induced 
by confined phonons as a function of energy splitting corre- 
sponding to the magnetic field region B £ [0.0, 2.36] T (before 
the anti-crossing) . Black solid/Green dashed curve represents 
the relaxation due to the deformation-potential/piezoelectric 
coupling. 



potential coupling and the piezoelectric coupling)^ the 
other one is the zeros due to the radial component, which 
induces different dips for different clcctron-phonon cou- 
pling mechanism. Note that the deformation-potential 
coupling dominates the relaxation rate in most cases ex- 
cept in the energy region [0.2,0.35] meV. This is because 
in this region, the corresponding phonon mode tends to 
be transverse mode which decouples to the electrons via 
the deformation-potential coupling.— 

It can be seen for Fig. [2] that the second feature, e.g., 
the local maximum of the relaxation rate, occurs in the 
region far away from the van Hove singularities, where 
the conhncd phonon DOS exhibits a plateau. In this re- 
gion, the behavior of the relaxation rate is decided by the 
combined effect of the ST mixing and the phonon emis- 
sion efficiency. Both mechanisms enhance the relaxation 
rate at the anti-crossing, resulting in the local maximum. 
For bulk phonons, the DOS is a quadratic function of the 
energy splitting, so the relaxation rate can be suppressed 
as energy splitting decreasing, resulting in the local min- 
imum at the anti-crossing. Note that the energy splitting 
E here (E = 0.14 meV) is much larger than the energy 
splitting in typical disk-shaped GaAs QDs, so the sup- 
pression of the relaxation rate due to the decreasing of 
the phonon emission efficiency is absent. 

To further justify the behavior of the relaxation rate 
at the anti-crossing, we give an estimation based on per- 
turbation theory. We concentrate on the deformation- 
potential coupling since it dominates the ST relaxation 



rate at the vicinity of the anti-crossing. We restrict the 
discussion with the lowest two basis loos') and |01T+) 
with energy eo an d ei + 5MB-B since they are almost de- 
generated at the vicinity of the anti-crossing. By ap- 
plying the degenerated perturbation theory, we get the 
eigenstates 



|S*) = cos6>|OOS)+sin0|OlT + ), 
|f+) = -sin0|OOS)+cos0|OlT+) 



(10) 
(11) 



where tan0 = t/(d + \J d 2 + 1 2 ) with t = 

7«.<#;(*)|9«l#(*)) and d = (e - ei + gfiaB)/2. 

The energy splitting between the two eigenstates is 
E = 2\ / d 2 + 1 2 . Note that the two states are maximally 
mixed as E reaches its minimum. 

The relaxation between the singlet state \S) and triplet 
state 1 7+) can be calculated by applying Eq. [9] For bulk 
phonons, we have 



r^ lk = e\mP 



ST 4 



bulk 



V 



Sill 



(I 



Iqt+ (Qc COS ( 



x e 



-(5/47r)(g c d o sin0) 2 



(12) 



^buikl = y/nD 2 /{Ad 3 z Tftv\p) is a constant. 



where 

q c = E/(Hvl) is the bulk phonon emission wave vec 



tor. 



W?) 



(00S\t 



+ e 



! |00S) - (OlT+le* 



e lqZ2 \0!T + ) is the axial component of the form factor. In 
the situation we considered here, 7g T (q) varies slowly, 
so we approximate it as a constant. For small E, one has 



]— ex E, 



(13) 



For confined phonons, only the lowest confined phonon 
mode has contribution at the vicinity of the anti-crossing, 
so we have 



Tconf rp— 1 I 7i rD I z ± 

ST+ — & iW conf l 



-^ST + (Qc 



l)q 



x li nt (q c , V2d Q )li nl (q c , do/x/2) 



<7=9c 
2 



(14) 



where |M° nf | = D 2 uj 4 /(2^/TTd z pR) is a constant. Note 
that the eigenencrgy of the first confined phonon mode 
tends to be linear in q, so we can write q c = E/(hvs)^ 
J| T (q) is the axial form factor which is the same as 
the bulk phonons since both the bulk and confined 
phonons can be expressed as plain wave in axial direction. 
I? onl (q,d ) = X $(k L /dQ) 2 e A ^y^J^~ z is the ra- 
dial form factor for confined phonons, with k^ = u vq jvL 

and A L (v,q) = (q 2 — fc|)d§/(47r). Xvq is the coefficient 
in the expression for the confined phonon cigenmode 
which is calculated numerically. It can be shown nu- 
merically that I^ onf (q,do) tends to a constant for small 
qM So for small E, we have 



r^conf 771- 

1 ST+ <X E 



(15) 
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FIG. 3: (Color online) Same as Fig. but for the ST_ relax- 
ation. 
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Comparing Eq. (|T3|) and Eq. (fl5|) . one can see that as 
the energy splitting E decreasing, the relaxation rate in- 
duced by bulk phonons decreases while the relaxation 
rate induced by confined phonons increases. 

Similar features also exists for ST_ and ST relax- 
ation rates. The multi-peak structure in the relaxation 
rate due to van Hove singularities can be found for ST_ 
relaxation rate shown in Fig. [3l Note that the "sharp" 
peaks here (the first, second and the fourth peaks) are 
not divergent peaks since they correspond to the q = 
van Hove singularities. They are "sharp" since the piezo- 
electric coupling has a large contribution at these peaks 
which can be seen in Fig. [HJd). This is because these 
peaks are axial mode as q — > 0, which is favored by the 
piezoelectric coupling^ Also note that for bulk phonons, 
the decreasing of the phonon emission efficiency due to 
the small phonon wave length dominates in this region. 
So although the DOS increasing as the magnetic field, 
the relaxation rate is suppressed. 

In the region where the confined phonon DOS exhibits 
plateaus, the effect of phonon emission efficiency can be 
pronounced, which may suppress the relaxation rate. As 
an example, we compare the STo relaxation rate induced 
by bulk and confined phonons in Fig. 01 For large mag- 
netic field, the relaxation rate is suppressed due to the 
decreasing of the phonon emission efficiency for confined 
phonons. For bulk phonons, the relaxation rate is still in- 
creasing due to the increasing of the phonon DOS. Note 
that the increasing of the relaxation rate for both con- 
fined phonons and bulk phonons in small magnetic field 
is due to the ST mixing. 



IV. SUMMARY 

In summary, we have investigate the ST relaxation in- 
duced by the confined phonons in two-electron nanowire- 
based QDs. We find that the behavior of the relaxation 
rate is dominated by the large DOS at the vicinity of the 
van Hove singularities, while in the region where the con- 



fined phonon DOS is flat, the effect of the phonon emis- 
sion efficiency and the ST mixing are more pronounced. 
This results in the multi-peak structure of the ST re- 
laxation rate and the local maximum of ST + relaxation 
rate at the vicinity of the ST + anti-crossing. These fea- 
tures are very different from the previous results for disk- 
shaped QDs, indicating the unique property of the con- 
fined phonons. The effect of the deformation-potential 
coupling and piezoelectric coupling are also discussed. 
We find that for InAs [111] nanowire, the deformation- 
potential coupling dominates the relaxation in most cases 
except in the region where the longitude component of 
the confined phonons mode is suppressed. The piezoelec- 
tric coupling is also found to be important for axial con- 
fined phonon mode at q — > 0. It is also worth noting that 
the relaxation rate induced by confined phonons is much 
smaller than the one induced by bulk phonons in most 
cases, this suggests the nanowire-based QDs are prefer- 
able for the application in quantum information and com- 
putation. 
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Appendix A: W? q /P (r) in Eq. [9] 



D /P 

The quantity W uq (r) is the radial component of the 

D/P 

clcctron-phonon coupling H C p with the superscript "D" 
for the deformation-potential coupling and "P" for the 



piezoelectric coupling. They can be obtained by substi- 
tuting Eq. © into Eq. ©/Eq. Q. 

For the deformation-potential coupling, W®(r) is 
given by 



W»(r) = (k 2 L +q 2 )J (k L r) X (? q \ 



(Al) 



^ 2 -q 2 . 



where k 2 L/T = {u uq /v L/T ) 2 

For the piezoelectric coupling for wurtzite InAs [111] 
nanowire, W{, q (r) is given by 

W* g (r) = [r'dr'Ioiqr^Koiqr^ 
,(o) 



x g 



xi U ^o(^r')(2g 2 -3fci) 



+ X^Mk T r')(k 2 T - Aq 2 )/q 2 ] /V3. (A2) 



where r< = min(r, r') and r> = max(r, r'). 
Jo(r)(Io(r) / K (r)) is the zeroth-order (Modified) Bcsscl 

function. xi>q an d Xi"j are coefficients in the expression 
for the confined phonon eigenmode which are calculated 
numerically following Refs. I38H33 . 
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